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Abstract. 

We present azimuthal anisotropy measurements in Au+Au and Cu+Cu collisions 
^ [ at ysjvjv = 62.4 and 200 GeV. Comparison between reaction plane and cumulant V2 

■ measurements in Au+Au collisions at y / 3jvjv = 200 GeV show that non-flow contributions, 

originating mainly from jets, influence the extracted V2 for pr > 3.5 GeV/c. Number of 
constituent quark (NCQ) scaling of V2, when studied as a function of transverse kinetic 
energy KEt, is seen to hold for Au+Au collisions at ^/sjvjv = 62.4 and 200 GeV and for 
Cu+Cu collisions at y / ijvjv = 200 GeV for KEt ;$ 1 GeV/c. Differential hexadecupole 
flow t»4 seems to exhibit scaling with integral t>2 for centrality < 40% as has been observed 
for differential V2- 
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1. Introduction 

Experimental measurements at the Relativistic Heavy Ion Collider (RHIC) have indicated 
the creation of a deconfined phase of quarks and gluons pQ. The azimuthal anisotropy of 
particles emitted suggests that the high-energy density matter rapidly thermalizes and its 
evolution is driven by large pressure gradients (2j [3] . Azimuthal correlation measurements 
in Au+Au collisions at RHIC have been shown to consist of a mixture of jet and harmonic 
contributions The azimuthal anisotropy extracted from these correlations are thus the 
sum of contributions from collective flow and jet-like correlations. The latter constitute one 
source of non-flow correlations, especially at high pr, and one would like to study how they 
vary with px and centrality. It has been reported that constituent quark scaling of V2 as a 
function of KEt or mj- — m is observed in Au+Au collisions at ^snn= 200 GeV and that 
differential V2 scales with integral i>2 for both Au+Au and Cu+Cu collisions at ^/snn= 
200 GeV |5] . These observations have implications on the hydrodynamic description of the 
medium as well as the pertinent degrees of freedom in the flowing matter j2j . PHENIX has 
measured the vi of charged pions, kaons and protons in Au+Au collisions at -^/stvtv = 62.4 
GeV and Cu+Cu collisions at y / i7v r /v= 200 GeV. It is interesting to study whether NCQ 
scaling holds for these two colliding systems and if the scaling of differential vi with integral 
V2 also works for differential U4. 



2. Data analysis 

The results presented here use data obtained by the PHENIX collaboration for Au+Au 
collisions at ,/sJm= 62.4 and 200 GeV and Cu+Cu collisions at ^snn= 200 GeV. The 
azimuthal anisotropy relevant to this work were measured using two different methods: the 
reaction plane [61 [7] and the cumulant [8] method. The reaction plane method measures 
the correlations between the azimuthal angle of the particles detected in the central arms of 
the PHENIX detector and the azimuth of the reaction plane obtained with the two Beam- 
Beam Counters (BBC) located at | rj |~ 3 — 3.9 [9]. The cumulant method extracts V2 from 
multiparticle correlations of particles detected in the PHENIX Central arms [TU] . The large 
rj separation between the BBC's and the central arms of PHENIX means that the extracted 
anisotropy is much less affected by possible non-flow contributions, especially those due to 
jets [9l [TT]. The differential second-order cumulant V2 is determined with respect to the 
reconstructed v% integrated over px- Monte Carlo simulations have been carried out to 
ensure that the azimuthal anisotropy is reconstructed faithfully with the PHENIX detector 
acceptance [El [10] . The particles used to reconstruct the differential flow are excluded from 
those used in the integral flow reconstruction to avoid autocorrelations. The second order 
cumulant ^2{2} is extracted from two-particle correlations and thus includes correlations 
from particles originating from jet fragmentation. Hence, the second order cumulant V2 is 
influenced by non-flow correlations from jets as discussed in [8]. The results presented in 
the next section show how these jet correlations affect the azimuthal anisotropy in Au+Au 
collisions at tJsnn= 200 GeV. 



3. Results 

3.1. Comparison between V2 using different methods of measurement 
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Figure 1. Comparison between V2 from cumulant and reaction plane (RP) methods for 
different centralities in Au+Au at y/s^N = 200 GeV. The closed and open circles, labelled 
V2{RP} and V2{2}, are the cumulant and reaction plane results respectively. 



Fig. [T] shows the differential vi in Au+Au at ^Js~nn = 200 GeV as a function of px 
from the reaction plane and the cumulant methods for several centralities. The measured 



V2 from the two methods agree well up to px ~ 3.5 GeV/c but differ for higher pr, with 
the cumulant V2 being higher than the reaction plane V2- This could be attributed to the 
influence of jets on the extracted V2 which is more prominent in the cumulant method. One 
also notes that the same trend is observed at all centralities, even though the agreement 
between the two sets of results holds up to higher px for the 0-10% centrality bin. 

3.2. Number of constituent quark scaling 

Number of constituent quark scaling of i>2 in Au+Au at s/snn = 200 GeV works well when 
studied as a function of transverse kinetic energy, as shown in the left panel of Fig. [2] [5j. 
This scaling scheme has been extended to the t>2 of positive and negative pions, kaons and 
protons measured in Au+Au collisions at y / SftF/v = 62.4 GeV and Cu+Cu collisions at ^/snn 
= 200 GeV. The middle panel of Fig. [2] shows that the constituent quark scaling of v 2 works 
well for tt^, K and (p)p in Au+Au collisions at ^/snn = 62 A GeV. The same good scaling 
of V2 is observed for the smaller Cu+Cu system at y/sJw = 200 GeV, as shown in the right 
panel of Fig. [2j These observations suggest that partonic degrees of freedom exist in the 
matter formed at lower beam energy and in the smaller colliding system. 
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Figure 2. V2/n q as a function of KEx/n q for charged pions, kaons and protons in minimum 
bias Au+Au collisions at ^/snn = 200 GeV (left panel), in Au+Au collisions at ^/snn = 
62.4 GeV for centrality 10 - 40% (middle panel) and Cu+Cu at ^/sTW = 200 GeV for 
centrality - 50% (right panel). 



3.3. V4 of unidentified particles 

PHENIX has recorded high statistics Au+Au collisions during the 7th year of data-taking, 
enabling the detailed study of higher harmonics of azimuthal anisotropy as a function of 
px and centrality. It has been reported that V2 scales with the integral flow, which has 
been shown to be proportional to the eccentricity j5]. The use of the integral flow scaling of 
differential V2 has been used to study the thermalization of the matter at RHIC [5]. Since, 
in general, v n ~ v^ 2 , which yields t>4 ~ v\ for n = 4, differential would be expected to 
scale with the square of the differential V2 and the scaling of differential V4 with integral V2 
could be investigated. The left panel of Fig. [3] shows the pt dependence of V4 of unidentified 
particles for three centrality bins in the range 10-40 % and the right panel shows V4, scaled 



by the square of the p^-integrated vi for the corresponding centrality (labelled e in the 
figure), as a function of pr- The latter plot shows that V4 scales with the square of integral 
V2 for centrality < 40 % and px < 2.0 GeV/c. By plotting fa/e ) 2 on the same plot (open 
circles) one observes a good agreement between differential t>4 and vf, when they are both 
scaled by the square of integral V2- 
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Figure 3. Left panel: Differential V4 of unidentified hadrons as a function of pt (left 

panel) for different centralities. Right panel: v^/e as a function of px, where e represents 
V2 integrated over pr. For the sake of comparison, the open circles represent differential V2 
squared scaled by the square of the integral V2- 



4. Summary 

In summary, we have presented results of azimuthal anisotropy from different methods 
at RHIC. They show that non-flow correlations from jets may influence the measured 
anisotropy for pt > 3.5 GeV/c. Constituent quark scaling is observed to hold in Au+Au at 
yJSNN = 62.4 and 200 GeV and Cu+Cu at ^Jsnn = 200 GeV. This suggests that partonic 
degrees of freedom are manifest at the lower RHIC beam energy and in the smaller Cu+Cu 
colliding system. The differential hexadecupole flow in Au+Au at y/s^N = 200 GeV is 
observed to scale with the integrated V2 for centrality < 40% and pt < 2.0 GeV/c. 
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